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The QCD factorization theorem for diffractive processes in DIS is 
■ used to derive formulae for the leading twist contribution to the nuclear 

JL . shadowing of parton distributions in the low thickness limit (due to the 

(— | \ coherent projectile (photon) interactions with two nucleons). Based on 



X 



the current analyzes of diffraction at HERA we find that the average 
strength of the interactions which govern diffraction in the gluon sector 
at x < 10~ 3 ,Qo = 2GeV is ~ 50mb. This is three times larger than 
in the quark sector and suggests that applicability of DGLAP approxi- 
mation requires significantly larger Qq in the gluon sector. We use this 
information on diffraction to estimate the higher order shadowing terms 
due to the photon interactions with N > 3 nucleons which are important 
for the scattering of heavy nuclei and to calculate nuclear shadowing and 
Q 2 dependence of gluon densities. For the heavy nuclei the amount of the 



*On leave of absence from PSU. 
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gluon shadowing: Ga{x,Q 2 ) )/AGn(x,Q 2 ) )\ x<1q s ~ 0.25 — 0.4 is sensitive 
to the probability of the small size configurations within wave function of 
the gluon " partonometer" at the Qq scale. At this scale for A ~ 200 the 
nonperturbative contribution to the gluon density is reduced by a factor of 
4 — 5 at x < 10~ 3 unmasking PQCD physics in the gluon distribution of 
heavy nuclei. We point out that the shadowing of this magnitude would 
strongly modify the first stage of the heavy ion collisions at the LHC en- 
ergies, and also would lead to large color opacity effects in eA collisions at 
HERA energies. In particular, the leading twist contribution to the cross 
section of the coherent J/ip production off A > 12 nuclei at y/s > 70 GeV 
is strongly reduced as compared to the naive color transparency expecta- 
tions. The Gribov black body limit for F2a(x, Q 2 ) is extended to the case 
of the gluon distributions in nuclei and shown to be relevant for the HERA 
kinematics of eA collisions. Properties of the final states are also briefly 
discussed. 
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I. INTRODUCTION 



It has been realized by Gribov already before the advent of QCD that there exists a 
deep relation between the phenomenon of high-energy diffraction and the nuclear shad- 
owing phenomenon M. In particular, the nuclear shadowing due to the interaction of 
a virtual photon with two nucleons can be unambiguously calculated in terms of the 
7* + N — > X + N diffractive cross section if the coherence length l c = Q2+° M 2 is much 
larger than the nucleus radius, Ra- Here M 2 is the invariant mass squared of the quark- 
gluon system to which a virtual photon is transformed [TJ. In the case of the charged 
parton structure functions (F 2 a{x, Q 2 )) connection between shadowing and diffraction 
has been explored for a long time, see [f2|-|TT|l and references therein. The importance 
of the color fluctuations- weakly interacting configurations in the shadowing phenomenon 
was first understood in H where this effect has been estimated based on the QCD aligned 
jet model and included in the calculation of F 2 a(x, Q 2 ). 

Additional contributions to the nuclear structure function are related to the piece of 
the photon wave function for which the coherence length l c is of the order of the av- 
erage internucleon distance r^N ~ 1.1 Fm. These important nuclear effects have been 
estimated and explored in P,[5|,|T2"|| using constraints which follow from the QCD momen- 
tum and baryon sum rules. Account for these effects leads to a more complicated QCD 
evolution which mixes shadowing region and the region of larger x. 



In the recent paper we started analysis of the implications of the information 
which is now available from HERA on the role of the gluon degrees of freedom in the 
diffractive processes in DIS for the gluon nuclear shadowing. We were able to study shad- 
owing for x < 1CT 3 and Q 2 ~ 20 — 50GeV 2 and predict a factor ~ 2 — 3 larger shadowing 
for the gluon channel than for the quark channel. This is in line with expectations of 



14], though it differs from the pattern assumed in a number of the models, see e.g. for 



the recent summary fT5|j . 

In this paper we will extend this analysis to a broad range of x and Q 2 . The main 



tool we will use is the QCD factorization theorem for the hard diffractive scattering [16 
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see also |T7[]. Application of the QCD factorization theorem makes it possible to establish 
correspondence between the deuteron(nucleus) infinite momentum frame (IMF) and the 
rest frame descriptions and therefore to explore advantages of both descriptions. 

An evident advantage of the IMF description is the simple interpretation of the mo- 
mentum and the baryon sum rules. On the other hand space-time development of high 
energy processes and nuclear shadowing phenomenon have a more clear interpretation 



within the nucleus rest frame approach. Using the QCD factorization theorem [I6| and 
the Gribov analysis of nuclear shadowing we will derive the model independent expres- 
sions for the leading twist nuclear shadowing of parton densities in the case of coherent 
interactions with two nucleons (section 2). [] However the shadowing due to the inter- 
action with two nucleons cannot diminish the parton density by more than a factor of 
0.75 without introducing ghosts into the theory (see discussion in the end of section 3). 



Hence in section 3 we use the recent analysis ]n| of the HERA diffractive data to extract 
the information on the S-channel dynamics of diffraction which is necessary to calculate 
the effects of coherent interactions with N > 3 nucleons. We find out that this anal- 
ysis implies that the average strength of the interaction responsible for the diffraction 
in the gluon channel at the resolution scale Qo ~ 2GeV and x < 1CT 3 is very large: 
a e ff ~ 50 — 60mfe. Large value of the interaction strength could be related to the large 
cross section of the small color octet dipole interaction with a nucleon which is given by 

37T 2 A 
a ™color octetdi P ole",N(.Einc) = a s{Q )xG N (x,Q = — ), (1) 

where x = 2m < ^E — ■ This is a factor of 9/4 larger than for the case of "color triplet 
dipole" ]I7|]. If we take for A the value we estimated before for the color triplet case: 
X(x ~ 10~ 3 ) ~ 9, we find that the cross section is close to the ^-channel unitarity limit 
for the range of applicability of DGLAP approximation corresponding to Oinei > ®ei for 
x ~ 10~ 4 and Q 2 ~ 10 GeV 2 . 



For an early discussion of the general arguments for the presence of the nuclear shadowing in 
the leading twist and references see Q. 
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For the scattering of a system with a radius r <tir N the slope of the elastic scattering 
is given by the square of the two-gluon form factor F 2 g (t) rs exp(B 2g t) with B 2g ~ 
4GeV~ 2 . For this situation condition <Ji ne \ = o e \ corresponds to the effective cross section 
of 87iB 2g = 40mb. This value is close to the one which emerges from the analysis of 
the diffractive data where the size of diffractive system is smaller though not negligible 
as compared to the nucleon size. Note also that the smallness of the shrinkage of the 
diffractive cone for the J/ip elastic photoproduction (AB < !GeV~ 2 for y/s between 5 
and 200 GeV as compared to AB rs 3GeV~ 2 expected in the soft regime) indicates that 
perturbative physics occupies most of the rapidity range for Q 2 > 4GeV 2 for HERA 
energy range. 

It is worth emphasizing that inapplicability of the DGLAP evolution equation and 
possible closeness to the unitarity limit we discuss here are due to the growth of 
xGn(x, Q 2 ) generated predominantly by the logQ 2 terms in the DGLAP evolution equa- 
tions rather than solely by the ln(l/a;) terms which would be the BFKL approximation. 
Hence the pattern discussed here is qualitatively different from the BFKL scenario of 
high-energy dynamics. Indeed, large values of xG N (x,Q 2 ) (~ 10 — 20 in the small x 
HERA kinematics and growing with Q as ~ y/Q) emerge not because of long ladders in 
rapidity - the ladders contain no more than 2-3 gluons in the multi-Regge kinematics, 
but rather due to a large number of emitters at the lower resolution scale. Possible close- 
ness to the unitarity limit makes it likely that for moderate Q 2 < lOGeV 2 corrections 
to the DGLAP predictions for the nuclear shadowing would be rather large. This would 
primarily affect our predictions for moderate Q 2 since the information about the gluon 
induced diffraction is obtained predominantly at larger Q 2 and extrapolated to lower Q 2 
via the DGLAP equations. 

In section 4 we first analyze the dynamics of the fluctuations of the interaction 
strength (color coherence - color opacity and color transparency phenomena) and ex- 
plain that significant fluctuations of the strength of interaction should be present in 
particular due to the QCD evolution. Next we study the nuclear shadowing originating 
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from the interactions with N > 3 nucleons. We point out that the eikonal type ap- 
proximation seems reasonable in the soft QCD regime when the projectile wave function 
contains a large number of constituents. On the contrary in the PQCD regime where 
photon wave function is given by a qq dipole not more than two inelastic collisions are 
allowed by energy conservation law. Otherwise the energy released in the inelastic colli- 
sions calculated through the cuts of exchanged parton ladders will be larger than the sum 
of the energies of the colliding particles. Evaluation of a larger number of rescatterings 
in PQCD is beyond the scope of the naive semiclassical approximation and requires an 
accurate account of the space-time evolution of the scattering process, in particular a 
calculation of the NLO approximation to the photon wave function. We demonstrate 
that the N > 3 interactions are sensitive to the existence of the fluctuations of the in- 
teraction strength. This sensitivity is rather small for A ~ 12. For such A we predict 
significantly larger shadowing for gluons: Ga(x, Qq) / AG n (x , Qo)|x<io- 3 Q 2 =4GeV 2 

~ 0.7 

than for quarks: F 2 a(x, QD/AF 2 n(x, <2o)|a!<io-3,Q§=4GeV 2 ~ °- 85 - For lar g er A sensitivity 
to fluctuations steadily increases. However we find that the average interaction strength 
in the gluon channel is large at the normalization scale of Qo = 2 GeV so a significant 
nuclear shadowing of average and larger than average interaction strengths into the cross 
section is determined by the geometry of collisions and rather insensitive to the structure 
of the distribution over the strengths. As a result of shadowing of strongly interacting 
(nonperturbative ?) configurations, the relative contribution of the interactions with 
small a is strongly enhanced in the parton distributions in heavy nuclei. We estimate 
possible effects of the weakly interacting configurations and find that they may contribute 
up to 50 % to G J 4~2oo(a ; < 10~ 3 ,Qo ~ 2GeV). Q At the same time the fraction of the 
cross section due to weakly interacting configurations should diminish with decrease of 
x. 



2 We are indebted to A.Mueller who stressed the effect of filtering of the PQCD physics in the 
parton distributions in nuclei. 
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We want to stress here that the use of information on the diffraction in DIS at HERA 
allows us to take into account the nonperturbative effects in the gluon nuclear parton 
densities at the boundary of the QCD evolution. In the previous studies the gluon 
shadowing either was treated purely perturbatively as for example in the IMF model of 
McLerran and Venugopalan or it was introduced in a phenomenological way assuming 
similarity between the shadowing in the gluon and quark channels, see e.g. [f|,|l5|. Overall 
a currently popular scenario which is used in the discussion of the heavy ion collisions 
assumes that reduction of gluon densities is a relatively small correction, for the recent 
review and references see |20| . 

Next, we introduce the constraints on the gluon densities which follow from the mo- 
mentum sum rule and imply presence of the gluon enhancement at x ~ 0.1. Combing this 
effect with the quark and gluon shadowing for small x we calculate the x, Q 2 dependence 
of the leading twist nuclear densities. In the end of the section we also consider nuclear 
structure functions in the limit when the nucleus thickness is large enough so the black 
disk approximation is applicable. 

Obviously, the predicted large gluon shadowing has many implications for the vari- 
ous high-energy processes of scattering off nuclei. In section 5 we calculate the impact 
parameter dependence of the gluon shadowing and briefly analyze two phenomena: the 
emergence of the color opacity in the coherent production of J/ip and T-mesons from 
nuclei in the HERA kinematics, and the suppression on the minijet production in AA 
collisions at the LHC energies. We find both the color opacity effect and minijet suppres- 
sion to be very large. For example, for the lead-lead collisions we predict a suppression 
of the minijet production at p t = 2(3) GeV/c by a factor > 7(> 4). 

In section 6 we briefly discuss properties of final states and predict a dip in the ratio 
of the spectra of leading hadrons produced in the current fragmentation region in eA 
and in eN collisions for rapidities shifted from the maximum rapidity by In (jVl^j^j /fi 2 
where is the average diffractive mass 2 produced in eN scattering. 

In section 7 we compare our approach with some of the recent studies of the nuclear 
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shadowing. 



II. THE QCD FACTORIZATION THEOREM AND THE LEADING TWIST 
SHADOWING FOR THE PARTON DENSITIES. 

The studies of the diffraction production in hard processes lead to the introduction 
of the diffractive parton densities f® B (/3,Q 2 ,xip,t) with j3 = which represent the 
number densities of partons in the initial hadron, but conditional on the detection of the 
diffracted outgoing hadron B in the target fragmentation region with light-cone fraction 
1 — xp and fixed momentum transfer t. For example in the case of the diffractive process 
e + p — > e + p + X the diffractive structure function F 2 D which is introduced via 

^U-w^^-y^-^), (2 ) 

can be written as 

F^iP, Q 2 , zip, t) = ^2 e lPfa/ P (P, Q 2 , x ip, t) + HT corrections. (3) 

a 

In the case of the proton production this structure of the hard diffractive processes was 
first suggested in the framework of the Ingelman-Schlein model [2I|. Recently it was 
demonstrated [16] that the QCD factorization theorem is valid for the x, Q 2 evolution 
of these parton densities at fixed xjp, t. The evolution is governed by the same DGLAP 
equations as for the inclusive processes. The HERA data on diffraction in DIS are 
consistent with the dominance of the leading twist contribution except near the edge of 
the phase space (see discussion below). 

For the processes dominated by the vacuum channel the Gribov theory JT| unam- 
biguously relates diffractive processes in the scattering of a projectile off a single nucleon 
to the process of nuclear shadowing due to the interaction of the projectile with two 
nucleons. The simplest way to visualize this connection for example in the case of the 
scattering off the deuteron is to consider j*d scattering in the deuteron rest frame in the 
kinematics where l c ^> Rd (Rd is the radius of the deuteron). Due to the difference of the 
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spatial scales characterizing the deuteron and the soft QCD strong interactions, the dom- 
inant contribution is given by the diagrams where the photon dissociates into a hadron 
component before deuteron and then this component interacts with both nucleons. Let 



us use the AGK theorem |22| and consider the cut of the double scattering diagram cor- 
responding to the diffractive final state (Fig.l). This corresponds to the scattering off 
one nucleon in the \in) state and off the second nucleon in the {out\ state. The final state 
interaction between nucleons is accounted for as usual within the closure approximation. 
The interference between two diagrams results from the Fermi motion of nucleons in the 
deuteron since the spectator nucleon in the \in) state has to have a momentum equal to 
the momentum of the diffracted nucleon in the (out\ state. The screening effect is ex- 
pressed ultimately through —Ref 2 where / is the diffractive amplitude of the interaction 
of the probe with the nucleon as compared to \f\ 2 in the case of diffractive scattering 
off the nucleon. The real part of the diffractive amplitude is rather small and can be 
calculated from the information about the imaginary part of the amplitude. Thus the 
difference between \f\ 2 and —Ref 2 is small and easy to deal with. 

Hence we can apply the Gribov results for the scattering off the deuteron and nuclei to 
evaluate the shadowing contribution to the deuteron parton density of flavor j in terms 
of the corresponding nucleon diffractive densities (we consider only the Pomeron type 
contribution, so we do not distinguish diffraction of protons and neutrons) 

f jPH (x, Q 2 ) = f j/p (x, Q 2 ) + f j/n (x, Q 2 ) -V^J dxjpdtS(At)f? /N (/3, Q 2 , x Fl t) . (4) 

Here S(t) is the electromagnetic form factor of the deuteron, and —t = {k 2 + 
(x F m N ) 2 )/(l -xp), and 77 = (1 - (ReA di} / ImA dif ) 2 ) /I + (ReA dif / ImA di} ) 2 ). 

Similarly, in the approximation when only scattering off two nucleons in the nucleus 
is taken into account one can similarly deduce the expression for the shadowing term in 
terms of the parton densities 

j p roo roo rXQ 

f j/A (x,Q 2 )/A = f j/N (x,Q 2 ) - - / d 2 b / dz x \ dz 2 / dx P ■ 

■if IN (A Q 2 > x Pi t) \ k 2 =0 PaQ>, zi)PA(b, z 2 ) cos(x F m N (z 1 - z 2 )). (5) 



Here pA{ r ) is the nucleoli density in the nucleus normalized according to the equation 
/ pA{r)d 3 r = A. For simplicity we gave the expression for the limit when the slope of the 
dependence of diffractive amplitude on the momentum transferred to target nucleon, t, 
is much smaller than the one due to the nucleus form factor so that impact parameters 
of two nucleons are equal. Note that eq.|5] is similar to the corresponding expression for 
the shadowing in the vector dominance model, see eq.(5.4) in [^] since the space-time 
evolution of the interaction is the same in both cases. This leads to the same structure 
of the nuclear block, provided one substitutes the VDM expression for the longitudinal 
momentum transfer, q z = My/2v by the Bjorken limit value: q z = XjprriN- 

The crucial feature of eqs^,^ is that the parton densities which enter in the shadowing 
term evolve according to the leading twist evolution equations. When they are folded 
with a function of xjp which does not depend on Q 2 they retain this property. Since the 
QCD evolution of real and imaginary parts of hard amplitude is governed by the same 
evolution equation at sufficiently small x we investigate in the paper the fact that real 
parts enter into diffraction and into shadowing in a different way does not influence the 
QCD evolution. This proves that eqs.^^ correspond to the leading twist contribution to 
the nuclear parton densities. In the limit of the low nuclear densities egs. [||| provide a 
complete description of the leading twist nuclear shadowing. 

Obviously the derived equations could not provide a complete picture of the deviations 
of nuclear parton densities from the sum of the nucleon densities for all x. This is because 
the derived equations take into account the contributions related to the distances l c ^> Ra 
but not the ones related to the configurations with much smaller coherence lengths. 
The simplest way to estimate the corresponding additional piece is to apply the energy- 
momentum and baryon sum rules which are exact in QCD for the leading twist parton 
densities. Therefore to satisfy these sum rules the shadowing should be accompanied by 
an enhancement of some parton densities at higher x. This enhancement term has to be 
added to eqs. £|,§. If we introduce this term at a scale Ql for x > xo it would contribute for 
large Q 2 for much smaller x. Hence the Gribov type approximation becomes inapplicable 
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for fixed x and Q 2 — > oo 0. Below, to deal with the enhancement effects we will adopt 
the procedure of 0] in which these effects are estimated at a low normalization point and 
the subsequent evolution is dealt with by solving the DGLAP evolution equations. 

The range of the validity of approximation where interactions with N > 3 nucleons 
are neglected strongly depends on the strength of the corresponding diffraction channel. 
Hence in the next two sections we review the results of the recent analysis of the HERA 
diffractive data and build approximation for treating interactions with several nucleons. 



III. DIFFRACTION AT HERA AND SHADOWING IN THE LOW NUCLEAR 

THICKNESS LIMIT 

A. Gap probability for the gluon induced hard diffraction 

First, let us briefly summarize the results of the studies of diffraction in D1S which 
were performed over the last few years at electron-proton and proton-antiproton colliders 
and recast them in the form necessary for the studies of the nuclear shadowing phenom- 
ena. The data obtained at HERA include studies of the diffractive structure functions 
in 7* + p — > X + p scattering, production of dijets and charm in 7* + p scattering. The 
leading twist contribution appears to describe the data well except very close to the edge 
of the phase space where higher twist effects are important. So the factorization theorem 
for these processes |l6j seems to hold for the studied Q 2 range, see |1| for the recent 



analysis. 

In the practical applications an assumption is usually made that the semiinclusive 
parton densities at small values of xjp can be written as a product of a function of xp 
and a parton density which depends on (3 = x/xp and Q 2 . 



In principle one should also take into account the effects of nonnucleonic degrees of freedom 
(the large x EMC effect) but for any practical purposes this effect is negligible. 
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Hence for the sake of brevity we will refer to these densities as the parton densities in 
the "Pomeron". ^ As a result one can define a diffractive parton density at given x as a 
convolution of the so called Pomeron flux factor, f p /ip(xip) and corresponding Pomeron 
parton density for example for gluons: 

xgdif(x,Q)= f P /A x ip) x / x ipfg/ip( x / x iP,Q )dx F , (6) 

J X 

where x max is the maximal value of xp for which diffractive picture still holds. 

The important finding of the HERA diffractive studies is that fg/p(/3) ^> f q /p(/3) for 
a wide range of (3 (Similar trend was observed in pp collisions, see review in ||26||). For 
example, in the best global fit of the HERA diffractive data |18| (fit D): 



PfM*> Qo) = (9-7 ± 1.7)0(1 - /?), Y/ g/ T (3 ^f w °' 13 - (7) 

Let us consider probability of diffractive events where the proton remains intact for 
the hard leading twist processes coupled solely to the gluons. It can be defined as 

xg N (x,Q 2 ) 



4 Note that in difference from the usual parton densities which are process independent the 
"Pomeron" parton densities may depend on the target, on the mass of diffractively produced 
system etc. In particular, for small masses M 2 <C Q 2 contributions of the higher twist to 
diffraction become important which are proportional to xGn(x, Q 2 ) 2 and hence lead to intercept 



of the effective "Pomeron" ajp(O) > 1.20 [24]. In the analysis of |T8| this kinematics was 
excluded from the fit. At the same time for M 2 3> Q 2 intercept should be more close to 
Qijp(O) = 1.08 familiar from the soft QCD interactions. These are particular illustrations of 
the deep difference between the QCD factorization theorem and the Regge pole factorization 



[25]. Note also that the energy dependence of diffraction is different from that for soft hadronic 
processes. This is not surprising since the coherence length for the soft hadronic processes is 
significantly larger than that for soft hadronic processes in DIS. 
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Since this definition includes only the leading twist contribution into diffraction it effec- 
tively excludes the contribution of small masses to the diffraction which could originate 
from the higher twist effects. Since at small x ~ 10~ 3 and Q 2 ~ few GeV 2 the ratio of 
the quark and gluon densities in a nucleon is ~ |, and f g /jp(P) ^> f q /ip(f3) one obviously 
expects 



p !if(x,Q 2 ) = f± Pl f (x,Q 2 ) » P q dtf (x,Q 2 ). (9) 



We can quantitatively estimate Pj^x, Q 2 ) using the fit D of |Tj|. The analysis of |T8| 



chooses the initial conditions for the DGLAP evolution at Ql = AGeV 2 - see eq.^. We 
also take x max = 0.02 which is the highest xp for which the one has enough sensitivity to 
the gluon density in the "Pomeron". However for x <C x max the diffraction probability 
practically does not depend on the choice of x max . We find that 

Pdifi 1 ^ 4 <x<3- 10~ 3 , Ql) « 0.34 • (1 ± 0.15), (10) 

which is much larger than Pg ap (x,Qo) ~ 0.12. Total probability of rapidity gap 
which includes double diffractive events (proton dissociation) is larger by a factor 
~ 1.4. This factor can be estimated assuming the Regge factorization for t = : 
^dJ^^^j^+pyif I t=o ~ O-^ independent of the diffraction state X\, and taking into 
account that the slope of the t dependence in the double dissociation should be smaller 
by about a factor of two due to almost complete disappearance of the proton form factor 
in the proton vertex. However the cross section given by the HERA groups includes a 
small contribution of the proton dissociation of about 15% |27|. So effectively the scaling 
factor is smaller ~ 1.25. 

Thus in the gluon channel the ratio of total diffraction to total cross section reaches 
the value ~ 0.4 for Q = 2 GeV. Thus we conclude the ratio of single diffraction to total 
cross section in the gluon channel is close to that for pp collision (for the soft hadronic 
processes analogous quantity is the ratio of sum of the elastic and the single diffraction 
cross sections to the total cross sections). 



13 



The QCD evolution leads to a decrease of this probability since at larger Q 2 many 
small x partons originate from "ancestors" at Qq with x > x max which cannot produce 
protons with small xp < x max . So at Q 2 = 25GeV 2 , Pj^ drops to about 0.2. [] 



B. Implications for the S-channel picture of hard diffraction 

The large probability of diffraction in the gluon channel comparable to that in soft 
hadron interactions indicates that configurations involved have large interaction cross 
section (here we effectively switch to the ^-channel language of description of diffraction 
281 , 12911 ). We can quantify this by using the optical theorem | t=0 = to introduce 



the strength of interaction a e ff as f\ • Q 

167Td<7 dl f /dtu-r, „ , , , 

a eff (x,Q 2 ) = '-^ = P! lf (x,Q 2 )16nB (11) 

a 

The effective cross section a e ff(x, Q 2 ) characterizes within the Gribov theory the diffrac- 



tive rescatterings of the produced quark-gluon system, cf. eq.|12[ The results of the 
calculation are presented in Fig. 2 for Q = 2GeV and x max = 0.02 for the quark and 
gluon channels and show that a e ff for the gluon channel is of the order 55 mb for small 
x and about 3 times larger than for the quark channel. Q Large value of cr e ff can be 



5 Note that experimental studies of the "Pomeron" gluon densities are performed either at 
large virtualities of Q > 5GeV or via scaling violation of f^. for Q > 2GeV. So they cannot 
directly measure the large value of P^j ■ 

6 Here and below we neglect < 5% corrections due to the real part of the amplitude since other 
uncertainties in the input are of the order 15 — 20%. 

7 For the sake of simplicity we parameterize the t dependence of diffractive cross section as 
da di f /dt = da di f jdt {t=0 exp Bt. 

determination of a e ff requires the knowledge of the t-dependence of the diffraction. Exper- 
imentally it was measured for the process 7* + p — » X + p only and for relatively large 
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interpreted as an indication that the interactions in the gluon channel is remaining strong 
up to much larger virtualities than in the quark channel. This matches rather naturally 
with the perturbative QCD pattern of a factor of 9/4 larger cross sections for color octet 



dipole -nucleon interaction than for the color triplet dipole-nucleon interaction [ OJ] , [|T 
see discussion in the introduction after eq.l. 

C. a e ff and shadowing for small nuclear densities 

Now we are in a position to rewrite the results of section 2 for the parton shadowing 
in the limit of small nuclear densities using the notion of cr e //. This would allow us as 
a next step to go beyond the two nucleon approximation for the shadowing effects. For 
x <C x m ax, and not too large A such that l c ^> R A the cos(xjpM N (zi — z 2 )) factor in 
eq.^ can be substituted by one. Hence in this limit the amount of shadowing is directly 
proportional to the differential probability of diffraction at t = 0: 

xG A (x, Q 2 ) = AG N (x, Q 2 ) - 4tt f d 2 bT 2 A (b) J*™* \ t=0 f g (x/x F , Q 2 )d X]P , (12) 

where 

T A {b) = / dzp A {z,b), (13) 



(xjp) ~ 0.01. For this kinematics the fit of [18] describes the t-dependence well. The fit also 



assumes the rate of the diffractive cone shrinkage oc exp(2o' ln(l/xjp)) with a' from the soft 
processes. Due to a larger value of diffraction and hence a larger value of o e ff f° r the gluon 



channel the assumption of [18] that the slope for single diffraction in the gluon channel is at 
least as large as for the quark channel seems also very natural. One should however remem- 
ber that the lack of direct measurements of the i-dependence of the gluon induced diffraction 
introduces an additional uncertainty in the results of calculations. Overall our guess for the 
uncertainty in the value of the parameter a e ff for 10 -4 < x < 10 -3 is about 20%. 
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is the usual nuclear thickness function and Pa{^) is the nuclear density. For larger xp one 
has to take into account the damping factor due to the cos(x]pMn{zi — z^)) factor which 
originates from the nuclear form factor due to the longitudinal momentum transfer in the 
transition to the diffractive mass equal to m^xp. However with our choice of x max = 0.02 
this effect is small even for A ~ 200. Eq.|l2| leads to the shadowing proportional to <7 e ff'- 

1 - xG A (x, Q 2 )/AxG N (x, Q 2 ) = ^J- J d 2 bT 2 (b) (14) 

and hence predicts a factor ~ 3 larger shadowing for the gluon channel than for the quark 



channel for Q = 2GeV . This is in line with expectations of though it differs from 
the pattern assumed in a number of the models, see e.g. for the recent summary ||15|| . 

Note also that if one does not separate leading and higher twist contributions in 
the diffraction off a nucleon one can still use the experimental data about the total 
cross section of diffraction off a nucleon to calculate using eq. [12] the total amount of 
shadowing in the corresponding channel for the scattering off the deuteron (and nuclei 
in the approximation when the interactions with N > 3 nucleons are neglected). The 



simplest way to see this is to apply the the AGK cutting rules [0] which are valid for the 



scattering off nuclei. In particular, if the the higher twist effects due to interactions with 



two nucleons described by the Mueller and Qiu model |3T| were important in the nuclear 
shadowing at the normalization point they should be manifested as well in the diffraction 
off a nucleon. So, as far as the diffraction is described by the leading twist factorization 



approximation, eq. [12] leads to the DGLAP evolution of the nuclear shadowing. 

Note also that in the approximation when a probe (photon) may interact not more 
than with two nucleons there exists a relation between the shadowing for the total cross 
section and the partial cross section of inelastic processes with the multiplicity similar to 
the one in the inelastic ep scattering - a±: 

&tot — Vimp ~ ^double, °~1 = & imp ~ 4(7 'double-, (15) 

where Oi mp is the impulse approximation cross section and Odoubie is the screening cross 



section due to the interaction with two nucleons pz| . One can see from eq.|15] that 
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shadowing due to the photon interaction with two nucleons can diminish the total cross 
section by not more than a factor of 0.75 as compared to the impulse approximation 
without introducing ghosts into the theory: for A e ff/A < 0.75 the partial cross section <j\ 
would become negative |^2[.0 This implies that for A e ff/A < 0.75 shadowing interactions 



with a larger number of nucleons could not be ignored. 

IV. FLUCTUATIONS OF THE INTERACTION STRENGTH AND 
MULTINUCLEON SHADOWING 

A. Modeling effects of cross section fluctuations 



Due to the large value of cx e ff(x,Q 2 ) for the gluon channel, deviations from eq.|T2 
due to interactions with N > 3 nucleons become large already for A ~ 10. To account 
for these effects we address the Q 2 dependence of <r e ff(x,Q 2 ). Within the DGLAP 
approximation it basically reflects an influx to small x of configurations which at a lower 
resolution Q' correspond to configurations with larger x = x parent and hence with smaller 
<7 eff( x 'parent, Q') ■ Configurations which interacted strongly at Q ~ Q interact strongly at 
large Q as well, but they contribute smaller and smaller fraction of the total cross section 
relevant for the nuclear shadowing phenomenon at a fixed x. This pattern is the same 
as in the QCD aligned jet model 0. Since the gluon shadowing strongly reduces gluon 
densities already Q ~ Qq the deviations from the DGLAP equations for Q ~ Q due 
higher order order terms in nuclear parton density originating from the average masses 
of the diffractively produced system should be significantly smaller than in the model 
of |3T| where shadowing at the starting scale is neglected. One can speculate of course 
that these effects have already occurred between Q 2 ~ lGeV 2 and Q 2 = Qq. Higher 



9 T 



Note that in Ref. 33] where shadowing was calculated in the Mueller and Qiu model the 



values of A e jf/A as low as 0.5 were obtained. 
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twist effects are enhanced for the contribution of diffract ively produced system with the 
masses M 2 <C Qq. This interesting question is beyond of the scope of this paper. 

It is straightforward to take into account effects of the longitudinal momentum trans- 
fer in the diffraction [[4|,|T3||. However as we demonstrated in that these effects are 
important only for xjp > 0.03 for A ~ 200 and even for larger xjp ~ 0.05 for light nu- 
clei. Since we have chosen x max = 0.02 we can safely neglect this effect in the following 
discussion. In this approximation to account for the fluctuation effects it is convenient 
to introduce the probability distribution over the strength of interaction in the gluon 



channel - P g (cr). a e ff is expressed in terms of P g (cr) as |34] 



a eff = J daa 2 P g (a)/ J daaP g (a). (16) 



We obtain in the generalized eikonal approximation: 

G A (x,Ql) _ Jd 2 bdaP g (a)(2-2exp(-T(b)a/2)) 



(17) 



G N (x,Q 2 ) AJdaaP g (a) 
where T(b) was defined in eq.|T3|. Fluctuations lead to a decrease of shadowing effect 
as compared to the quasieikonal approximation where P{cr) oc 5(cr — er e //)H We will 
study effects of fluctuations at length elsewhere. However for characteristic a ~ a e ff ~ 



55 mb the exponential factor in the numerator of eq.[I7] is very small for A ~ 200 and 
small enough b a wide range of a, leading to cross section m 2-nR 2 A . This suppresses 
the contribution of large a -nonperturbative QCD physics and therefore enhances the 
contribution of small a-PQCD physics. 

As a result of large absorption for a ~ a e f / fluctuations near the average value of a 
practically do not change the shadowing (for A < 250) provided cr e // is kept fixed. If, for 
example, we assume that P(<j) = a9(cr — cr ) the Ga/Gn would change even for heavy 
nuclei by less than 20%. 



10 Note that since the value of <r e // is fixed by the cross section of the diffractive scattering 
this approximation differs from the eikonal approximation often used for the hadron-hadron 
scattering in which only elastic rescatterings are included. 
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Besides, as we mentioned in the introduction the diffraction in the gluon induced hard 
processes without proton break up constitutes about 35% of the events at x < 10~ 3 . Large 
contribution of small a in / daaP g (a) would imply that for scattering in a significant 
fraction of configurations the gap probability is much smaller than average. Hence some 
configurations would have to generate the gap events with a probability exceeding 35%. 
This is one of indications for the problems for the applicability of DGLAP to describe ep 
scattering at HERA at low normalization point Qq. At the same time there is a natural 
mechanism in PQCD for the generation of a contribution of very small a in / daaP g (a) . It 
comes from the QCD evolution. The partons at given x which originated from ancestors 
at lower resolution scale with Xi n iu a i > x m ax do no contribute to diffraction and hence 
effectively correspond to a ~ (In other words this effect reflects the contribution of 
small coherence lengths to the small x physics due to the QCD evolution). This effect is 
especially pronounced in the models of QCD evolution with low normalization point like 
the GRV model |35|| . To estimate the sensitivity to this effect we will use two models. 
The first one is the quasieikonal which neglects higher order fluctuations. The second 
model is the fluctuation two-component model (to which we will refer to as a fluctuation 
model) which implements an extreme assumption that a fraction A of the total cross 
section at given x originates from configurations with small cross section and the rest 
from the average ones. The second model is similar to the QCD aligned jet model of 
]2||§]. It corresponds to 

aP g (a)(xX5(a) + (l-X)5(a-a ), (18) 



where do is fixed by eq.16 to 



"° = (19) 

The requirement that the gap probability for this component (including dissociation of the 
nucleon) does not exceed 50% puts an upper limit on A. Taking into account uncertainties 
in the value of the total rapidity gap probability we estimate that A < 0.2. So we will use 
A = 0.2 in this model. Note that for small values of a e ff < 20mb and moderate values 
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of A the screening weakly depends on A. Hence we will ignore this effect for the quark 
channel. For simplicity we will also assume that A does not depend on x. This should be 
considered as a rather rough approximation since for small enough x < 0.005 one may 
expect the contribution of small l c to decrease with decrease of x. Also for larger x the 
drop of the parameter a e ff can be due to an increase of A. However in this region we 
anyway have significant uncertainties due to the contribution of xjp ~ x max . Thus the 
results in this x range can be considered simply as a smooth interpolation between the 
region of large shadowing and the region where shadowing disappears. 



Obviously, small cr's in eq.[T7] would give the dominant contribution for A — > oo. 
However for large cr e // and A < 0.2 for A ~ 240 a contribution of the weakly interacting 
component could at most become comparable to the soft contribution, see Figs. 3,4 below. 



B. Models for shadowing at Qq and numerical results for Q 2 dependence 

Based on the above discussion of the fluctuation effects we adopt the following pre- 
scription in the calculation for the gluon and quark nuclear shadowing for A > 10: 

(i) For Q = Q we use two models: the quasieikonal model and the fluctuation model 
with A = 0.2 analogous to the ones we used in The relation between the amount 

of shadowing in the two models is given by 

G A {x,Q )/G N (x,Q ) 

Qo) = 

A + (1 - X)G A (x, Q )/G N {x, Q ) quasieik.mod. ( T ' — ) (20) 

(ii) To constrain the behavior of the gluon density at x > 0.02 in the normalization point 
we use the analysis of || which indicates that gluons in nuclei carry approximately the 
same fraction of the momentum as in a free nucleon: 

/ dxxG A (x, Q 2 ) / dxxG N (x,Q 2 ). (21) 
Jo Jo 

This allows to estimate the amount of the gluon enhancement at x > x max assuming 
that it should be concentrated at x < 0.2 where average longitudinal distances (the 
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Ioffe distances) contributing to the parton density are comparable to the internucleon 
distances - this procedure is similar to one we introduced in J5|. 

(hi) Based on the rational presented above we use the DGLAP evolution equations 
to calculate the nuclear shadowing for larger Q 2 . 

Since in this paper we are interested primarily in the behavior of the nuclear gluon 
and quark densities at x < 0.01 we are not sensitive to details of the enhancement pat- 
tern. So we do not try to introduce an A dependent shape for the enhancement and 
assume that the enhancement is present for 0.02 < x < 0.2 and can be approximated 
by G A (x,Q )/G N (x,Qo) = C(A)(x - 0.02) (0.2 - x). We also do not model small en- 
hancement for F 2 a(x,Q 2 ) and Va(x,Q 2 ) at x ~ 0.1. In the calculations we use the 
standard Fermi step fit to the nuclear densities: p(r) —Cj{\ +exp((— R + r)/b)), where 
R = 1.1/m • A 1//3 , 6 = 0.56/m. First, in Fig.3 we present a comparison of the gluon 
shadowing calculated in two models at Qq = 2GeV and x = 10 -3 , 10~ 4 . One can see 
that shadowing for A ~ 12 is already very significant, though the effect of fluctuations is 
still small. With a further increase of A the effect of fluctuations becomes larger and it 
reaches a factor of 1.4 for A~ 200. For these A in the fluctuation model weakly interact- 
ing configurations contribute approximately half of the cross section. Note however that 
we expect that in a more realistic model of fluctuations a relative contribution of these 
contributions is likely to decrease with decrease of x, see discussion in section IV.A. 
In Figs.4,5 we present results for the Q 2 dependence of shadowing for gluons and for 
quarks for A=12, 40, 100, 200 calculated in the quasieikonal model. One can see that 
the gluon shadowing is large already for A=12 and for heavy nuclei reaches the level of 
A e ff/A ~ 0.2. (Variations of the parameter a e ff within a factor of 1 ±0.2 allowed by the 
uncertainties in the value of the gluon diffractive density at t — lead to the similar vari- 
ations of G A /AG N at x < 10~ 3 and of [1 - G A (x,Q 2 )/AG N (x,Q 2 )} for 0.02 < x < 0.2.) 
Shadowing decreases with increase of Q 2 but remains large up to very large Q 2 . The 
shadowing in the charged parton channel is much smaller and rather weakly decreases 
with Q 2 . In fact, at x ~ 10~ 4 shadowing for F 2a {x,Q 2 ) first increases with increase 
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of Q 2 due to a larger shadowing in the gluon channel. Note also that eq.|2T] leads to 
a rather large enhancement of Ga/Gn for 0.05 < x < 0.2. However the enhancement 
is rather large already for A = 12. As a result a further growth of the enhancement 
between carbon and tin is of the order of 10% and it is well consistent with the analysis 
of the current data in |36 |. The decrease of shadowing with increase of Q 2 is due to the 
feeding of the small x by partons which originated from x > 0.02 at Qq. Effectively, as 
we discussed above, the QCD evolution leads to fluctuations in the value of interaction 
strength due to the mixing of the contributions of small and large l c . Therefore the use 
of eq.[l7| in the quasieikonal model with o oc <7 e // defined at a high resolution Q 2 would 
lead to an overestimate of the nuclear shadowing since the cross section fluctuations lead 
to a decrease of the shadowing for the fixed value of cr e ff-> see discussion in (T3I . 



Note also that it is often stated in the literature that nuclear shadowing for the 
total cross section of should be practically the same for the real photon and for 
virtual photons with moderate (few GeV 2 ) virtualities. However our analysis predicts 
a significant drop of the nuclear shadowing in the total cross sections of 7* A scattering 
between Q 2 = and Q 2 = AGeV 2 due to the strong Q 2 dependence of the diffraction 
contribution to a r *N- The results of our calculation using data on the diffraction in the 



7P scattering at W ~ lAGeV |37j and at HERA []38[ and a smooth interpolation between 
two energy ranges is presented in Fig. 6. For a rough estimate we use the quasieikonal 
approximation which leads to a slight overestimate of shadowing for the heavy nuclei. 



C. Structure functions of nuclei in the the black disk limit 

It is of interest to consider also the structure functions in the limit Q 2 =const, x <C 
m R , A — > 00 first analyzed by Gribov for ot-i&l 0- If wa s argued in this paper that in 
such a limit interactions of all essential configurations in the virtual photon wave function 
with a nucleus can be treated in the black body (^-channel unitarity) limit. Since in 
the black disc limit the dispersion over the strengths of interactions can be neglected one 
finds 
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gjM*. Q ) = — ^ (m2 - Q2)2 , (22) 

and 

1 2 vrQ 2 ^ /-as^f p{m 2 )dm 2 

& FlaM) = ^2^L (m 2 + Q 2 ) 2 ' (23) 

where p(m 2 ) = <j(e + e~ — ► hadrons)/a(e + e~ — > fi + fi~). 

We can generalize the Gribov formulae deduced for the sea quark distribution [I] to 
the case of the gluon channel: 

*GU(s, Q )/Q = -^-j ^ p(m ) (m2 + Q2)2 (24) 

Here p(m 2 ) = cr(J — > hadron) / o(j — > 2g) is the ratio of the cross section of the gluon 
hadronic processes initiated by the local operator J = -j=F^ x F£ x introduced in to 
the perturbative cross section of the annihilation of the gluon source into gluons. For 
m 2 — > oo, p{m 2 ) = 1. So in the black disk limit 

xG A (x, Q W = IS M = In (20 , (25) 



3vr 2 y Q 2 2m N R A J 8vr 2 Vx 

where Xo = • Eq.|25] provides a solution for the problem of the gluon nuclear 

shadowing in the theoretical limit Q 2 =const, x <C xo, ^4 — > oo. It also illustrates that 
investigation of the Q 2 dependence of the nuclear structure functions will provide an 
effective method to establish whether fluctuations of strengths play a significant role in 
the gluon structure functions. 



It is useful to compare eq.[25| with the impulse approximation value of Ga = AG^. For 
example for A = 200 this equation leads to xG^lO -4 , 10GeV 2 )/A = 8 while the current 
fits to the nucleon data lead to xGat(10~ 4 , lOGe^ 2 ) ~ 20. Hence the model independent 
unitarity constrain implies a large gluon shadowing for this kinematics. This is consistent 
with our model calculations presented in the previous subsection. 



23 



V. ONSET OF COLOR OPACITY REGIME IN HARD DIFFRACTION AND 
SUPPRESSION OF MINUET PRODUCTION IN NUCLEUS-NUCLEUS 

COLLISIONS 



Let us now briefly discuss some of the consequences of the found magnitude of the 
gluon shadowing. 

The production of minijets is often considered as an effective mechanism of producing 
high densities in the head on heavy ion collisions. However in the LHC kinematics for 
the central rapidities minijets are produced due to collisions of partons with xj et = Jg^ N ■ 
For heavy ion collisions s^n > 4 TeV and the gluon-gluon collisions are responsible for 
production of most of the minijets. Therefore the gluon nuclear shadowing would lead 
to a reduction of the rate of the jet production due to the leading twist mechanism by 
a large factor up to pt ~ lOGeV/c, see Fig. 7 where we give results of calculation in the 
quasieikonal and fluctuation models. The nuclear gluon shadowing leads to a similar 
very strong reduction of the heavy onium production in pA and AA collisions at LHC 
energies for y c , m . ~ and small p t . 

For the central impact parameters the reduction is even larger. Using generalized 
eikonal approximation of eq.O we can calculate also the suppression of the parton den- 



sities at a given impact parameter b as |40 |: 



r ( n 2 h , I daP g {a){2 - 2 exp(-T(b)a/2)G N (x, gg) 

J r g \a)oaa 

Since in the impulse approximation GA(x,Q 2 ,b) imp = GN(x,Q 2 )TA(b) we finally obtain 



The results of calculation using eq.^7| are presented in Fig.8. 

So we conclude that the shadowing effects are likely to reduce very substantially the 
parton densities generated at the first stage of heavy ion collisions at the LHC energies. 
However theoretical uncertainties related to the role of point-like configurations lead to 
rather large uncertainties in the estimate of the suppression for the case of the heavy 
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ion collisions. Hence it would be very important to perform a direct measurement of the 
gluon nuclear shadowing at HERA in this kinematics. 

Presence of a large gluon shadowing leads to large effects in the diffractive eA collisions 
at HERA energies. Here we consider the simplest example - coherent diffractive produc- 
tion of vector mesons at large Q 2 by the longitudinally polarized photons - 72+^4 — > V+A, 
and photo(electro) production of heavy onium states. Since these processes are domi- 
nated by production of qq in a small size configuration one may naively expect that color 
transparency should hold for such processes and the amplitudes of these processes should 
be proportional to A. However the QCD factorization for exclusive processes leads to 



the amplitude of this process been proportional to Ga(x, Q ) pi] . For small x and t 

dt ~ a[) G%(x,Q 2 ) dt ' 1 ' 

So we expect that the color transparency regime for x > 0.02 (with a small enhancement 
at x ~ 0.1 due to enhancement of Ga(x, Q 2 ) for these x) would be followed by the 
color opacity regime for x < 0.01. As an illustration in Fig.9 we present the ratio of 
the cross section of the J/ip and T production off nuclei with A = 12,200 and nucleon 
calculated under the assumption that the leading twist gives the dominant contribution 
in these processes. It is plotted as a function of x and normalized to the value of the 
ratio at x = 0.02. In the calculation we use the analysis of |^T| which indicates that 
Qlff ~ 5(40)GeV 2 for J/ip(T) production. One can see from the figure that we expect 
onset of the Color Opacity regime for J/i[) starting at x ~ 0.01. The Color Opacity effect 
remains quite significant for production of T. 



VI. NUCLEAR EFFECTS IN THE INCLUSIVE LEADING HADRON 
SPECTRA IN EA COLLISIONS 

Numerous data on hadron-nucleus scattering at fixed target energies indicate that the 
multiplcities of the leading hadrons Na(z) = — ^ aA , da ^ + dz — — decrease with increase of 
A. Here z is the light-cone fraction of the projectile "a" momentum carried by the hadron 
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"h" . On the contrary, the QCD factorization theorem for the inclusive hadron production 
in DIS implies that in the case of electron-nucleus scattering no such dependence should 
be present in the DIS case. This indicates that there should be a interesting transition 
from the soft physics dominating in the interactions of real photons with nuclei to the 
hard physics in the inclusive hadron production in the DIS kinematics. It would be 
manifested in the disappearance of the A-dependence of the leading spectra at large z: 

N A (z, Q 2 ) = N N (z, Q 2 ), for z > 0.2, Q 2 > few GeV 2 , (29) 

At small x a new interesting phenomenon should emerge due to the presence of diffraction 
and nuclear shadowing for smaller z. Indeed, the diffraction originates from the presence 
in the wave function of 7* of partons with relatively small virtualities which screen the 
color of the leading parton(partons) with large virtuality and can rescatter elastically 
from a target (several target nucleons in the case of nuclear target). Inelastic interactions 
of these soft partons with several nucleons should lead to a plenty of new revealing 
phenomena in small x DIS eA scattering, which resemble hadron-nucleus scattering but 
with a shift in rapidity from y max (current) related to the average rapidities of these soft 
partons. This shift can be expressed through the average masses of the hadron states 
produced in the diffraction: 

Vsoft partons ~ Vmax ~ hl^M^ //J 2 ), (30) 

where \i ~ lGeV is the soft scale. Partons with these rapidities will interact in multiple 
collisions and loose their energy leading to a dip in the ratio r)A(y) = Na(v) /N p (y). At 
the same time these multiple interactions should generate a larger multiplcities at smaller 
rapidities. Application of the AGK rules indicates that for y < y so ft partons — A, where 
A = 2 — 3 the hadron multiplicity in the case of nuclei will be enhanced by the factor: 

AF 2p (x,Q 2 ) 

MV) = (31) 
At the rapidities close to the nuclear rapidities a further increase of r]A(y) is possible due 

to formation of hadrons inside the nucleus. A sketch of the expected rapidity dependence 

of r) A (y) is presented in Fig. 10. 
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One also expects a number of phenomena due to long range correlations in rapidity. 
This includes: (a) Local fluctuations of multiplicity in the central rapidity region, e.g. 
the observation of a broader distribution of the number of particles per unit rapidity, due 
to fluctuations of the number of wounded nucleons ||. 

These fluctuations should be larger for the hard processes induced by gluons, for 
example the direct photon production of two high p t dijets. (b) Correlation of the central 
multiplicity with the multiplicity of neutrons in the forward neutron detector, etc. 

Another important manifestation of nuclear shadowing is a large probability of diffrac- 
tive final states for small x. In the case of the generic e + A scattering this probability 
would reach ~ 35% for A ~ 200 ||. Due to a larger effective cross section of interaction 
in processes induced by the hard interactions with gluons, this effect should be even more 
pronounced in the lepton scattering processes where a hard process corresponds to hard 
7* — g interaction. Hence for example in the charm electroproduction of a heavy nucleus 
we expect about a half of the events to originate from the coherent diffraction where 



the nucleus remains intact [ 13 1 . Also one expects more strong filtering out of the gluon 



dominated diffraction than the diffraction dominated by the coupling to the quarks. 

VII. COMPARISON WITH OTHER APPROACHES 

Several approaches were developed over last decade to the dynamical calculation of 
nuclear shadowing phenomenon in DIS R 

First group of approaches is based on the Gribov work [[jj which established connection 
between the diffraction and nuclear shadowing. Among these considerations pH^PlJTl 



the one of JTTJ] is the most detailed and comes closest to our analysis in the case of F2A- 
The analysis is based on the fit to the HERA diffractive data on the e + p — > e + X + p 
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A phenomenological approach based on fitting the existing nuclear data and imposing the 



momentum and baryon sum rules in the spirit of || was pursued in [12,1 
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reaction within the model developed by the authors and a set of assumptions about 
higher order screening effects which are anyway rather small for i 7 ^. They demonstrate 
that the model can well describe the NMC data at x ~ 0.01 and give predictions for the 
HERA kinematics. The main differences from our approach are the use of the model for 
diffraction which does not explicitly satisfy the QCD factorization theorem for diffraction 
in DIS and neglect by the effects of enhancement of gluon distributions in nuclei at 
x ~ 0.1. Besides the gluons appear to play a rather small role in their model of ep 
diffraction leading to the expectation of the gluon shadowing smaller than in the case of 
F2A as compared to the larger shadowing for gluons expected in our analysis. 

In the case of hadron-nucleus scattering both total cross section and inelastic diffrac- 
tion can be described based on the idea of the fluctuations of the interaction strength in 
the projectile treating interaction of each component in the eikonal approximation, for 
the review and references see |f£2" |. There are a several models where a similar approach 
has been applied to the calculation of the nuclear shadowing by introducing the impact 
parameter qq virtual photon wave function \l/ 7 (&) and introducing the cross section of 
the qq-N interaction for the fixed b, see e.g. |43| , fP| and references therein. So far it was 
assumed in these models that Uqq-N = cb 2 . hence the QCD evolution which leads to 
a fast increase of u q q-N with incident energy was neglected. Moreover in this approxi- 
mation shadowing for the small b configurations is a higher twist effect, leading to an 
expectation of lack of the leading twist shadowing for and lack of gluon shadowing. 

Another group of approaches uses the infinite momentum picture and treats all the 
process within the perturbative QCD. To avoid problems with positivity of the cross 
section for A e ff/A < 0.75 (see discussion in section XX) one has to include interactions 
with N > 3 nucleons. The current models which include such interactions (see |45j-47| 
and references therein) assume that all shadowing is generated perturbatively and do not 
include information about diffractive processes. Qualitative expectations of these models 
are a rather large gluon shadowing and significant nonlinear effects in the evolution of 
the parton densities. It would be interesting to compare two approaches after the lead- 
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ing twist shadowing effects are implemented and constraints following from the HERA 
diffractive data are taken into account. 

VIII. CONCLUSIONS 

We have demonstrated that dominance of gluons in the " Pomeron" diffractive parton 
densities leads to a large enhancement of the nuclear gluon shadowing in a wide range 
of x, Q. Gluon shadowing of this magnitude will strongly affect the first stage of the 
heavy ion collisions at LHC, lead to a number of Color Opacity phenomena in the HERA 
kinematics for the eA collisions. Study of the gluon shadowing for heavy nuclei may allow 
to enhance contribution of the small interaction strengths, allowing to unmask PQCD 
physics in the eA collisions at HERA. The studies of coherent diffraction off nuclei and 
hadron production in the inclusive eA scattering will provide complementary handles for 
studying the small x dynamics. 
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FIGURES 




FIG. 1. Double scattering diagram for the j*D scattering corresponding to a diffractive 
final state 
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FIG. 2. Dependence of c e // as denned in eq.[ll] for gluon and quark channels on x for 
Q = Q = 2 GeV. 



34 




35 




FIG. 4. Dependence of Ga/AGn on x for Q=2,5,10 GeV (dashed, dotted, solid curves) 
calculated in the quasieikonal model. 
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FIG. 5. Dependence of F2A/AF2N on x for Q=2,5,10 GeV (dashed, dotted and solid curves) 



calculated in the quasieikonal model. 
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FIG. 6. Dependence of A J ^ on for Q=0,2 GeV (dashed and solid curves). The data 



points are from [p9| . They are corrected for the small effect of the nuclear shadowing in the 
deuteron. 
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FIG. 7. Suppression of the jet production in AA collisions due to gluon shadowing at y = 
calculated in the quasieikonal and fluctuation models (solid and dashed curves). 
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FIG. 8. Impact parameter dependence of the gluon shadowing for the scattering of Pb for 
Q = 2GeV and x = 10~ 3 , 10 -4 calculated in the quasieikonal and fluctuation models. 
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FIG. 9. Color opacity effect for the ratio of the coherent production of J/tp and T from 
carbon(lead) and a nucleon normalized to the value of this ratio at x = 0.02 calculated in the 
quasieikonal and fluctuation models (solid and dashed curves). 
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FIG. 10. Sketch of the A dependence of the hadron multiplicity ^ for small x plotted as a 
function of hadron rapidity. 
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